Spin-glass behavior in Ni-doped Lai.gsSro.isCuC^ 



O 

(N 



ON 
(N 

O 

o 
u 

Oh- 
o3 



c 

o 
o 



> 

On 

NO 

in 
o 



X 



A. Malinowski, V. L. Bczusyy, R. Minikaycv, P. Dziawa, Y. Syryanyy, and M. Sawicki 
Institute of Physics, Polish Academy of Sciences, 02-668 Warsaw, Poland 
(Dated: January 18, 2013) 

The dynamic and static magnetic properties of Lai.85Sro.i5Cui_j,Ni a 04 with Ni concentration up 
to y=0.63 are reported. All the features that characterize the spin-glass (SG) behavior are found: 
bifurcation of the dc susceptibility \ vs temperature curve, a peak in the zero-field cooling branch 
of this curve accompanied by a step in the imaginary part of the ac susceptibility \ a c, frequency 
dependence of the peak position in the real part of Xac and scaling behavior. The decay of remnant 
magnetization is described by a stretched-exponential function. The characteristic time from the 
critical slowing-down formula that governs the dynamics of the system suggests the existence of 
the spin clusters. The strongest interactions between the fluctuating entities are at y equal to the 
charge carriers concentration in the system. The SG transition temperature decreases linearly with 
decreasing y for y<0.30 and extrapolates to K at y=0 what means that Ni display a magnetic 
character in the surrounding Cu-0 network starting from the smallest concentration y. The static 
critical exponents characterizing the scaling behavior of the nonlinear part of x he between those 
typical for three dimensional Ising-like and Heisenberg-like systems. 

PACS numbers: 74.72.-h, 74.62. Dh, 75.50.Lk, 75.40.Gb 



I. INTRODUCTION 

Transition-metal oxides (TMO) seem to be a na- 
ture's playground where a strong interplay between spin, 
charge, and lattice degrees of freedom appears. When 
doped with holes, their phase diagram involves spin- and 
charge-density-waves, Jahn- Teller distortions, ferroelec- 
tricity and superconductivity, to mention only few of 
their complex ground states^ The spin-glass (SG) phase 
is also quite common to the phase diagram of TMO 
and has been observed, among others, in cobaltites^ 
manganites,— £ nickelates^ and cuprates.— & This can be 
viewed as a natural consequence of a disorder present in 
such complex systems due to their chemical doping with 
heterovalent cationsJ^— 

In the phase diagram of cuprates, the region display- 
ing the SG features separates the antiferromagnetic (AF) 
and superconducting (SC) regions and its characteristics 
are observed even inside the SC dome^ ' 13 ' 14 Glassines of 
spins, accompanied by charge glass-like behavior, seems 
to emerge with the first added holesi 15 ' 16 At low temper- 
atures, these holes localize, presumably at oxygen sites, 
and form a local singlet on the Cu sites The mag- 
netic frustration may be understood as a result of lo- 
cal ferromagnetic exchange coupling between Cu 2+ ions 
in the background of three-dimensional (3D) long-range 
AF order^ The specific role of holes, including their 
self-organization in the stripes, and the details of the AF 
order, including the possible existence of the spin spi- 
rals, are still under debatei 11 ' 19 ' 20 The early observations 
of the spin dynamics and its evolution with doping have 
indicated that the entities undergoing the SG transition 
are the finite-size AF domains, cooperatively freezing in 
the inhomogeneous but magnetically ordered phase^^i 
Thus this phase is often refereed to as a cluster SG. While 
the thermodynamical SG-type behavior here is well es- 
tablished, the microscopic origin is still discussed 



The essential nature of cuprates is believed to be gov- 
erned by the physics of the Cu-0 planes, which are 
sometimes regarded as ideal 2D Hciscnberg spin systems. 
However, any practical realization of such idealized lay- 
ers is connected with introducing some inevitable, inter- 
nal, disorder! 12 ' 23 Although the sources of such a disor- 
der can be cuprate-family specific (such as, for example, 
the tilted CmOq octaedra and the stripe structure around 
£=0.12 in La2- a; Ba a: Cu04), the random distribution of 
Coulomb impurity potentials in the layers separating the 
Cu-0 planes is the most common origin. Introduction of 
heterovalent dopants or extra oxygen atoms is necessary 
to change and control the charge-carriers concentration 
in cuprates but creates a modification of the Coulomb 
potential that disrupts the lattice periodicity and thus 
serves as an additional scattering center for carriers in 
the Cu-0 planes^ 

It is natural that the above-mentioned internal 
Coulomb-driven disorder is minimal in these compounds 
in which the carrier reservoirs are located far from the 
Cu-0 layers. For example, the compounds such as 
HgBa 2 Cu0 4 and YBa 2 Cu 3 6+ y (YBCO) are close to 
the cuprate clean limit i 10 ' 23 ' 24 Doping with Ca allows to 
move the YBCO apart from the clean limit in a con- 
trolled way. The recent muon-spin-rotation (/iSR) ex- 
periments on this system give the arguments against the 
picture of the cluster SG phase introduced by disorder 
and suggest the common ground state, named frozen 
antiferromagnet (FAF), both for low-doping regime (i.e. 
in the re-entrant AF phase) and for larger dopings 
for which clusters of spins coexist with percolating 
superconductivit y 22 ' 23 The study shows that the local- 
field distribution is narrow, in the strong contrast to the 
canonical SG in which the distribution width is compa- 
rable to the mean value. The disorder does not mod- 
ify significantly the hole-concentration dependence of the 
FAF transition temperature. The origin of FAF remains 
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unknown but dilution of magnetic moments, not frustra- 
tion, seems to play the crucial role^ 

The situation in the dirty cuprates, such as 
La 2 _ 2; Sr 2 ;Cu04 (LSCO), is even more complicated be- 
cause of the intrinsic disorder. In LSCO, in contrast to 
the clean- limit compounds, the cluster SG "order" gov- 
erns solely the physics of the system in a quite large 
region of the phase diagram. The hole-doping depen- 
dence of irreversibility in magnetization, which is found 
in LSCO, has been interpreted in terms of a generic quan- 
tum glass transition , 15 ! 25 ' 26 it has also been suggested to 
be linked to the presence of the SC correlations above T c ^- 
and, recently, has been re-interpreted as an argument 
towards existence of spin-density-wave quantum critical 
point i2I The NMR measurements have revealed the en- 
hancement of spin- freezing temperature, T g , deeply in 
the SC phase of LSCO at xc^.0.12 and insensitiveness of 
this phenomenon to 1% of atomic disorder^ Whether 
the physics underlying the cluster SG phase in the non- 
superconducting region of LSCO phase diagram and that 
at x~0.12 is the same, is an open problem. The observa- 
tion of charge glass-like behavior in lightly doped LSCO 
rises a question of how this dynamic charge order evolves 
into SC state with increasing xJ£ Clearly, the interplay 
between superconductivity and magnetic and electronic 
glassiness in LSCO is far from being well understood. 

Impurities intentionally introduced into the Cu-0 
planes have been widely employed to probe the prop- 
erties of cuprates, with the hope that the response of the 
system reveals the generic features of the pure compound. 
So far, such a situation close to being perfect has been 
found only in the ID correlated systems. The magni- 
tude of magnetization induced in quasi- ID AF spin chain 
Y 2 BaNi05 by different nonmagnetic impurities appeared 
to decay exponentially, with the correlation length £j mp 
equal to the numerical prediction for the spin-spin £ in 
the pure systemJ^ In quasi-2D cuprates, the spinless de- 
fects in the Cu-0 planes induce paramagnetic moments 
on the surroundings Cu ions and disorder driven by this 
local magnetism extents the SG region in the phase dia- 
gram in comparison to the pure system. As a result, the 
phase diagram of Zn-doped clean-limit YBCO becomes 
similar to that of non-doped dirty-limit LSCO^ This 
may be viewed as an argument that intrinsic or extrinsic 
- magnetically driven in this case - disorder influence the 
system in the similar wayJ^ 

In the context of the above, Ni is an unique dopant 
in cuprates. This results from its ability to built NiOg 
octaedra - increasing the substitution of Ni ions into the 
Cu-sites in LSCO leads eventually to isostructural nick- 
elate La2- £! ;Sr a ;Ni04 (LSNO). Although superconductiv- 
ity is not observed in LSNO, its phase diagram is still 
very reach. The static ID charge ordering in the form 
of stripes, predicted in the context of high-Tc supercon- 
ductors (HTSC)22n2i has been first observed in LSNO^ 2 . 
and only later in Nd-doped LSCO^ The proximity of 
the charge ordering region to the SG phase suggests that 
they may be relatedi^i Thus, a systematic examination 



of the evolution of the LSCO properties with the Ni dop- 
ing may shed some light onto the problem of possible 
correlation between superconductivity, stripe order and 
the SG phase in HTSC. 

A remarkable evolution of views on the role of Ni 
dopant in cuprates has taken place, owing to accumula- 
tion of new experimental facts. Although nominally mag- 
netic (3d 8 , S=l), Ni 2+ ion appears to have weaker effect 
on superconductivity than the nominally non-magnetic 
Zn 2+ (3d 10 , S=0) ion^-il The dc susceptibility, x, mea- 
surements have revealed very small (0.6 /ib) paramag- 
netic moment of Ni introduced into the Cu-0 planes^ 
It has been suggested that at small concentration Ni 
is substituted as Ni 3+ ion^ The normal-state electri- 
cal transport experiments have suggested that the quasi- 
particle scattering at the Ni impurity has predominantly 
non-magnetic character.— On the other hand, the c- 
axis optical conductivity measurements in undcrdoped 
NdBa2Cu30g.8 show strong enhancement of the normal- 
state pseudogap energy by the Ni doping while doping 
with Zn is found to suppress the pseudogap. This dif- 
ferent impact on pseudogap has been attributed to the 
magnetic character of the Ni dopant^ 

Just recently, it has been claimed that the Ni does 
not disturb the AF spin-1/2 network in the Cu-0 planes 
when its concentration, y, is smaller than the hole con- 
centration, p, in the system42 The careful measurements 
of the local distortions around the Ni ions replacing Cu 
ions in La2- a; Sr x Cui_j ) Niy04 suggest that for small con- 
centrations, y < p (where p is equal to Sr content, x), 
Ni serves only as a hole-absorber and creates a strongly 
hole-bond state, called Zhang-Rice doublet^ with the 
effective moment S=l/2 that couples AF with S=l/2 mo- 
ments of the surrounding Cu ions. These measurements 
have been carried out on the single crystals in a wide 
p-range (jp=x, 0<x<0.15) but within limited Ni concen- 
tration (j/<0.07). Based on them, a magnetic-impurity 
picture for Ni dopant in superconducting cuprates has 
been completely disqualified, at least below the optimally 
doped regions 42 

Our study shows that the actual situation is much 
more complicated. We have carried out the dynamic and 
static magnetic measurements of polycrystalline samples 
of Lai.85Sr .i5CuO4 (LSC015) doped with Ni up to large 
concentration ?/=0.63, exceeding p (equal to x=0.15 holes 
per C11O2 plane) over four times. We have found the 
SG behavior in all non-superconducting samples, even in 
these with small y, below p. The transition tempera- 
ture T g , when extrapolated inside the SC region, takes 
finite values and approaches zero in the y=0 limit. Thus, 
the magnetic role of Ni should not be neglected even for 
y < p. Despite the microscopic mechanism, the low- 
temperature phase of Ni-doped LSCO exhibits all fea- 
tures characteristic for the SG systems: irreversibility in 
the T-dependence of Xi a peak in the zero-field cooling 
(ZFC) branch of the x(T) curve accompanied by a step 
in the imaginary part of ac-susceptibility Xac, the moder- 
ate (typical for the cluster SG) frequency dependence of 
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the peak position in the real part of Xaa and the scaling 
behavior. 



II. EXPERIMENTAL DETAILS 

The policrystalline samples of Lai.ssSro.isCui-yNiyC^ 
(LSCNO) with 0<y<0.63 were synthesized from 4N-5N 
pure La2C>3, SrCC>3, CuO and NiO by using conventional 
solid-state reaction method. The stoichiometric amounts 
of the powders were carefully mixed, pressed into pellets 
and sintered in a pure oxygen gas flow at 1320 K for 48 
hours. After cooling down to room temperature with the 
rate 2 K/min, the samples were reground and the whole 
procedure was repeated two times. 

The X-ray diffraction measurements were carried out 
at the Bragg-Brcntano diffractomctcr (XPcrt Pro Alpha! 
MPD, Panalytical, with a setting described in Ref. |43|). 
For determination of absolute value of unit cell parame- 
ters, NIST SRM 676 Alumina was used as an internal ref- 
erence material. Crystallographic characterization and 
structure refinement was done with help of FullProf.2k 

44 

program^ 

The magnetic susceptibility measurements were car- 
ried out at two setups: the commercial SQUID magne- 
tometer (MPMS, Quantum Design) working in the tem- 
perature range 2 K - 400 K and field up to 5 T (dc and 
ac measurements) and the home-build setup based on 
Cryogenics Consultants SQUID sensor and working in 
the temperature 1.6 K - 280 K and magnetic field up to 
0.3 T (a part of the dc measurements). 



III. RESULTS 
A. Crystallographic analysis 

Structural and phase analysis has shown that all sam- 
ples have tetragonal K2NiF4-type structure (space group 
I4/mmm) and include no impurity phases. An example 
of Rietvcld refinement plot is shown in Fig.[T]for LSCNO 
with the Ni content y=0!9. 

The lattice parameters of the undoped 
Lai.85Sro.i5Cu04 sample are a=3. 77709(8) A and 
c=13. 2368(4) A; thus, the axial ratio c/a is equal to 
3.5045(2). The powder diffraction patterns indicate 
absence of any structural phase transition in the whole 
investigated range of the Ni content. The lattice 
constants change linearly with increasing y: c decreases 
with the rate Ac/Ay = -5.9(1) x 10~ 3 A (wO.02% 
c) per 1 at.% of Ni while a increases with the rate 
da/dy = 8.1(1) x 10~ 4 A («0.05% a) per 1 at.% of 
Ni (see Fig. As a result, the c/a ratio decreases 
with the rate about 0.07% per 1 at.% of Ni. These 
results are consistent with the results of the earlier 
studies for a large Ni-doping range ; 45 ' 46 but show a 
considerably smaller scatter of the experimental points. 
The dc/dy and da/dy rates found here for y < 0.63 
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FIG. 1: The x-ray powder-diffraction pattern of 
La1.85Sro.15Cuo.8iNio.19O4. The observed pattern is de- 
noted by dots while the calculated one is marked as solid 
line overlaying them. The bottom curve shows the difference 
between them. The short vertical lines mark the positions of 
Bragg reflections. 
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FIG. 2: Lattice constants of LSCNO at room temperature as 
a function of Ni doping. 



are comparable to the rates found for y < 0.06 in the 
detailed X-ray-absorption fine-structure (XAFS) and 
X-ray powder diffraction (XPD) measurements reported 
in Ref. [47| . 

Free coordinates of the atoms in the tetragonal unit 
cell of LSCNO have been determined with the use of the 
Rietveld refinement procedure. Position (0,0, z) of the 
apical oxygen atom 0(2) as a function of Ni content is 
discussed in Appendix [A] Unlike the XAFS technique, 
XPD is not atom-specific in the sense that only aver- 
aged interatomic Cu/Ni-0(2) distance can be obtained. 
However, the detailed analysis (see Appendix[A) suggests 
that Ni-0(2) distance may change between y=0.07 and 

y=i- 

The bond length between the Ni atom and the in-plane 
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T(K) 

FIG. 3: The low-temperature dc-susceptibility of LSCNO 
measured at 10 Oe for various Ni content, y. As an example, 
the ZFC and FC branches of x(T) curve and the characteristic 
temperatures, T g and Ti rr , are indicated for y=0.63 sample 
by arrows. 



oxygen atom 0(1), Ni-O(l), has been found by XAFS 
measurements to have two distinct values, depending on 
the ratio of the hole concentration to the Ni content^ 
As mentioned, the XPD technique is not capable of re- 
vealing such atom-specific details. Since the 0(1) po- 
sition, (^,0,0), is constant in K 2 NiF4-like structure of 
LSCNO, the averaged value of Cu/Ni-0(1) distance, re- 
sulting from the XPD analysis, increases with y in pro- 
portion to lattice parameter a, i.e. with the same rate 
Aa/(aAy). This means contraction of Cu(Ni)Og octa- 
hedra along the c axis with increasing Ni content, ac- 
companying by simultaneous expansion in the ab plane. 
The resulting volume change in the whole doping range is 
very small, of the order of 0.1%. More details on analysis 
procedure, free atomic positions and bond distance will 
be given elsewhere 4& 



B. Normal-state dc susceptibility 

The measurements of x(T) at 10 Oe reveal that the 
superconductivity survives in LSCNO up to y=0.054. 
Above this Ni concentration, the normal-state x vs T 
curves display bifurcation below a characteristic temper- 
ature, Ti rr , depending on the thermal-magnetic history 
of the sample, i.e. whether it was cooled down in the zero 
magnetic field (ZFC mode) or in the non-zero field (FC 
mode). The representative curves are shown in Fig. [3] 
The ZFC branch has a rounded but very well defined 
peak, position of which, T g , is slightly below Ti rr - Let us 
note that in the canonical SG usually an opposite situa- 
tion is observed, i.e. Ti rr < T g . The behavior observed in 
LSCNO can be a manifestation of the existence of mag- 




T(K) 

FIG. 4: The dc-susceptibility of LSCNO, measured at 1 kOe 
field, for Ni content, y, equal to (from bottom up) 0.09, 0.11, 
0.13, 0.15, 0.17, 0.19, 0.25, 0.30, 0.50 and 0.63. The solid lines 
are the best fits to the Curie- Weiss law in the 30 K - 400 K 
range. Inset shows the characteristic temperatures Tc and T g 
and effective fitting parameter 9 (left scale) as a function of 
y. The solid line is the linear fit to T g vs y dependence for 
y < 0.30. In addition, the fitting parameter C is shown as 
a function of y in the whole doping range (right scale). The 
dashed lines are guides to the eye. 



netic clusters in the system.— 

The temperature of the peak in the ZFC x(T) curve, 
T g , which we take as the temperature of a transition to 
the SG phase, increases linearly with y up to y=0.30 
(see the inset to Fig. |4}. The smallest concentration y, 
for which we observe the bifurcation of x{T) curve, is 
equal to 0.056. The temperature T g , when extrapolated 
linearly outside our temperature measurement window 
and inside the SC phase, has finite values and decreases 
to zero at y=0 (see the solid line in the inset to Fig. [4}. 
This strongly suggests that Ni dopant in the sublattice of 
surrounding Cu spins displays a magnetic nature starting 
from the smallest concentration y. 

It should be emphasized that the deviation of x(T) 
measured in the FC mode from that measured in the 
ZFC mode cannot alone prove that a SG state develops 
below the bifurcation temperature of x{T) curve. The su- 
perparamagnetic systems and conventional ferromagnets 
with a broad distribution of potential barriers exhibit 
such a feature as well. In addition, the ZFC-FC splitting 
of x{T) is observed in disordered AF compounds that 
may be more relevant to the LSCNO system^ Thus, in 
the next sections, we will present the additional char- 
acteristics confirming the transition to the SG phase: a 
logarithmic frequency dependence of the real part of ac 
susceptibility, x! , below the transition, a step in the imag- 
inary part, x": a small frequency dependence of the peak 
position in x'(T) curve, described by the standard critical 
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slowing-down formula, time decay of the thermorcmnant 
magnetization and, finally and the most conclusive, the 
scaling behavior of the nonlinear dc-susccptibility. 

Before that, however, we will shortly describe the x be- 
havior at temperatures well above T g . The normal-state 
x(T) at larger temperatures for Ni concentration below 
y=0.09 can be analyzed in terms of the universal empiri- 
cal curve F, proposed for LSCO without parametrization 
independently by Johnston^iand Nakano.— The univer- 
sal function F was found to describe the effective suscep- 
tibility of the Cu 2+ spin sublatticc in LSCO with various 
Sr content when the x(T) data are represented in the 
reduced parameters (x - Xo)/(x ~ Xmax) and T/T max , 
where xo is the T- independent sum of three components: 

(1) - the isotropic contribution from the closed Cu shells, 

(2) - the Van Vleck contribution and (3) - the contribu- 
tion from the charge carriers, and Xmax is the maximal 
value that x(T) reaches at T max . Using the universal 
function from Ref. [39|, we have fitted x(T) to the for- 
mula x(T) = xo + A ■ F(T/T max ) + C/T and have cal- 
culated the effective magnetic moment introduced by Ni 
ion, fJ-eff, from the Curie term, C/T = Nfj/^^f/SksT, 
with the rest of parameters having their usual meaning 
(see Ref. 52]). The moment /U e // is constant up to 
y=0.07 and equal to (0.7 ± 0.05)/ie per Ni ion, what is 
very close to the value found previously when a linear 
function instead of F has been used in the analysis.— 
Note that there is no need to add any finite paramagnetic 
Curie temperature, 9, in the Curie term to reproduce the 
experimental data well, as it was necessary for Ni-doped 
LSCO with x=0.18 Sr content^ 

The x(T) term described by the F function is not 
detectable when Ni concentration exceeds y—0.09. For 
these large y, the x(T) data above ^30 K are well re- 
produced by the Curie- Weiss (CW) law, with the finite 
negative 9, x — Xo + C / {T — 9). The best fits are shown 
by the solid lines in Fig. |4j The calculated magnetic 
moment fJ> e ff does not remain constant, as it does for 
smaller y, and increases with increasing y to reach about 
1.6 fin per the Ni ion for y=0.5. Assuming temporary 
that all Ni ions are in the same magnetic state this would 
mean that the spin of Ni ion is much closer to 5=1/2 
(for which ^, e ff~1.73 jUg) than to 5=1 (/i e //=2.82/^). 
Note that the identical value (1.6±0.1)/Lts was obtained 
for Ni-doped YBa2Cu306+x, where Ni moment remains 
quasi-constant with hole doping in the investigated range 
up to y=0.04^ 

The CW term in the temperature dependence of sus- 
ceptibility indicates the presence of the localized mag- 
netic moments in the system. For small y<0.09, when 
the magnetism of the Cu-ions sublattice is clearly de- 
scribed by the function F and it is not much altered 
by the Ni ions introduced in small quantities, the Curie 
term in x(T) can be fully attributed to Ni ions, all be- 
ing in the same magnetic stated The resulting /i e // = 
gy / S(S + 1) = 0.7 fiB means that the spin 5 of Ni ion, 
immersed in the surrounding Cu ions sublattice, is equal 
to 0.11. Such a small value can be explained by the for- 



mation of the strongly hole-bound state of Ni ion, rep- 
resented as Ni 2+ L, as proposed in Ref. 42|. However, 
we observe the increase of the AW/ starting at j/«0.09, 
not at y=0.15 (being equal to the hole concentration in 
the system) , as it is expected from the picture suggested 
in Ref. |42|. Abrupt increase of the C in the CW law 
for t/>0.07 suggests that at least part of the Ni ions is 
in the magnetic state different from Ni 2+ L, even at con- 
centration below y=0.15. However, the macroscopic sus- 
ceptibility measurements alone do not allow to separate 
out the observed effective moment into two components 
coming from the Ni ions in two different magnetic states. 

The values of 9 for a given y display some dependence 
on the T-range used in fitting x{T). This makes to treat 
9 as the effective parameter rather than the true con- 
stant. In addition, there might be some crystal electric 
field effects, averaged in the policrystalline samples. Nev- 
ertheless, the negative sign of the effective paramagnetic 
CW temperature 9 indicates that the dominant exchange 
interactions in the system are antiferromagnetic. As it is 
shown in the inset to Fig. 21 the 9 exhibits strong de- 
pendence on Ni content. Its absolute value, \9\, increases 
abruptly when y increases from y«0.09 up to j/=0.15 This 
is followed by a rapid decrease, and a regime of satura- 
tion for y>0.30 (with a weak tendency to decrease with 
increasing y). This behavior can be understood as a re- 
sult of trapping the mobile holes by the Ni ions and will 
be discussed in the Sec. IIV CI 



C. Ac susceptibility 

The mctastable SG state is usually characterized by 
ac susceptibility, Xac- In Fig. [5] we show Xac{T) curves 
for LSCNO for several Ni contents, y=0.17 (LSCN017), 
y=0.25 (LSCN025) and y=0.50 (LSCNO50). A quite 
sharp cusp is visible in the in-phase component, x'{T), 
for all samples. Its temperature, Tf, is always a bit larger 
than the temperature of the maximum in dc~x(T), T g , 
and roughly coincidences with the temperature of the in- 
flection point in the step of the out-of-phase component, 
x"(T). At larger temperatures, above Tf, x"(T) is equal 
to zero, while below Tf has a finite value. Such a be- 
havior is characteristic for SG transition and allows to 
distinguish the SG compounds from the disordered AF 
systems, in which x"(T) is constant and remains equal to 
zero even below the temperature of the transitioni 50 ' 54 i 55 

In Fig|B] we depict the real part of Xac, normalized 
to the value at 1 Hz, as a function of frequency for 
LSCN025, as an example. At 25 K, i.e. above T g , in the 
paramagnetic state, the variation of / over three decades 
does not influence x' i n a noticeable way. Below T gi 
in the SG state, the x' exhibits a logarithmic frequency 
dependence. Such a frequency dependence has been pre- 
dicted theoretically for a short-range Ising SG^ and it 
has been observed for many SG systems^ However, this 
logarithmic relationship is not unique for the SG, because 
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FIG. 5: The real x' an d imaginary x" components of the 
ac susceptibility of LSCNO with y=0.17, 0.25 and 0.50 Ni 
content as a function of temperature, measured at f=l Hz 
and with 1 Oe amplitude of the applied ac field and without 
any dc-bias. 



it is exhibited by any two-level disordered system, dy- 
namics of which is governed by activated processes with 
a broad distribution of the activation barriers heights.— 
Note that at 8 K, closer to T g , the system is more sen- 
sitive to variation of frequency then at 2 K, "deeper" in 
the frozen state. 

The ac measurements at various frequencies reveal that 
the position of the maximum in x'{T) curve, 7/, moves 
toward higher temperature and the magnitude of x' de- 
creases with increasing frequency. Such a behavior is 
expected for a SG system. As a raw measure of this 
frequency dependence, the parameter S = dhiTf(uj)/duj 
is used^ In experimental practice, a shift between two 
outermost frequencies accessible in experiment is usually 
employed and S is calculated as S = ATj / (TfA\nuj)i 49 ^ 55 
The values of 5 obtained in this way for LSCNO system 
do not show any obvious correlation with y and are equal 
to 0.012, 0.012 and 0.014 for j/=0.17, 0.25 and 0.50, re- 
spectively. 
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FIG. 6: Normalized real part of Xac for LSCN025 as a func- 
tion of frequency for different temperatures - above and below 
T g . The amplitude of ac field was 1 Oe and the measurements 
were carried out in 10 Oe dc field. The error bars for the fre- 
quency scan at T=25 K are similar to these for the scans at 
2 K and 8 K. 



For systems with noninteracting entities (which can be 
cither particles or magnetic clusters), the values of 6 of 
the order of 0.1 have been reported^ and such a rela- 
tively high sensitivity to frequency is predicted in the 
classical model of superparamagnetism i 61 i 62 Any inter- 
actions between particles weakens this high sensitivity 
and this is why almost two orders of magnitude smaller 
values of S are measured in the canonical spin glasses 
(e.g. (5w0.002 for CuMn, Ref.[l|) or in the SG phase 
of other systems, such as manganites, where o«0.003 is 
reported.— Finally, in the case of well-ordered ferromag- 
netic or AF systems with even stronger interactions only 
MHz and GHz frequencies are sufficiently large to cause 
any observable shift in the x'(T) peak position.— 

The values of 5 ~ 0.01 obtained for LSCNO lay be- 
tween these extremes and are typical for cluster glasses, 
i.e. systems with randomly distributed interacting mag- 
netic clusters^ Recently, S = 0.016 was observed in 
Cc2CuGe3, which is an example of so called nonmag- 
netic atom-disorder SG system with possible (ferromag- 
netic) clusters^ To recapitulate, the use of the simple 
but model-independent criterion based on the parameter 
8 strongly suggests the possibility of the existence of spin 
clusters in LSCNO. The frequency dependence of Tf will 
be discussed in more detail in Sec. IIVBI 



D. Magnetic hysteresis 

The presence of spin clusters should be reflected in 
the magnetic hysteresis. Evolution of the magnetiza- 
tion loops in LSCNO50 is illustrated in Fig|7l where the 
isothermal M (H) curves are depicted at several temper- 
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FIG. 7: Evolution of the M vs H behavior with temperature 
for y=0.50 sample. The M — H isothermal curves are shown 
for 2 K (well below T g =13.1K), 10 K (still below T g ), 17 K 
(between T g and Ti rr ) and 50 K (well above Tf). Insets show 
details at low fields at T=10 K and at T=17 K. 



atures. The curves for LSCN017 and LSCN025 have 
qualitatively similar features. No trace of saturation is 
visible even at the lowest T. At 2 K, the individual 
branches of the M(H) curve obtained after ZFC have 
a characteristic " S" shape - the initial slope of the curve 
is smaller than the slope at inflection point at nonzero 
field. This is particularly clearly visible in the virgin 
curve. Such "S" shape of M vs H curve is typical for SG 
systems in a frozen stated With increasing temperature 
the S-shaped curve smoothly evolves into a straight line, 
indicating the paramagnetic behavior at high T . Hystere- 
sis is observed only below T g . The coercive field, taken 
from the hysteresis loops measured at 2 K, increases lin- 
early with y and is equal to 1.0, 1.4 and 2.6 kOe for 
y=0.17, 0.25 and 0.50, respectively. A high-field part of 
the virgin curve (above 9 kOe for y=0.50 at 2 K) lies 
outside the hysteresis loop. This suggests the presence of 
metastable states in the system. 

The existence of a hysteresis loop clearly excludes su- 
pcrparamagnctism as a candidate for the ground mag- 
netic state of LSCNO, since superparamagnetism is a 
thermal equilibrium behavior^ 7 - In the case of ferromag- 



netic or canted AF clusters existing in the system, a hys- 
teresis is expected but with large initial susceptibility be- 
cause clusters are at first saturated along their local easy 
axis and only after that various clusters became fully 
aligned along the applied field 42 If clusters reverse their 
magnetization coherently then hysteresis loops have al- 
most rectangular shapes^ In LSCNO just the opposite 
behavior is observed - the initial x is smaller then this 
at inflection point in M(H) curve. Thus, it is reason- 
able to conclude that M(H) curves do not show any fea- 
ture indicative of (large) clusters in LSCNO. The M(H) 
curves remain non-linear for T g < T < Ti rr indicating 
that SG phase starts to built over a broad T-range. For 
LSCN025, nonlinearity in the M vs H dependence is 
observed even at 25 K, i.e. at T ~ 2.5T S , well above 
Ti rr . The deviations of the M(T) curve measured at 
constant field from the CW dependence fitted in a wide 
T-range (up to 400 K), appear at around 30 K. These 
features indicate that short-range AF correlations and 
fluctuations exist well above T g . Such precursors of co- 
operative freezing at T g have been observed in metallic 
SG with 9>0 at temperatures higher than 5Tj£ and in 
the amorphous SG with both signs of 8 at even higher 
temperatures T«12 - 20T g M^ 



E. Decay of thermoremnant magnetization 

Since the SG system in a frozen state can react to the 
applied field slowly, the magnetization curves obtained 
after ZFC procedure do not give information about the 
thermal equilibrium of SG but rather reflect a slow in- 
crease of magnetization (and thus susceptibility) with 
time 42 The underlying physics is similar to the one that 
governs the decay of the remanent magnetization with 
time. 

To investigate the mechanism by which the system de- 
cays back to equilibrium we have applied a following pro- 
cedure. First, the magnetic field of 1000 Oe has been 
turned on at 200 K. Subsequently, the sample has been 
cooled down to 200 K during 80 min, then the tempera- 
ture has been kept constant for 10 min, and, finally, the 
magnetic field has been switched off. Next, the remnant 
magnetization has been measured vs time at T = 2 K, 
starting immediately after the field became zero. As it 
can be seen in Fig [51 the time decay of this FC ther- 
moremanent magnetization, Mtrm, is well described by 
a stretched exponential formula, 



TRM 



(i) = MoexpHt/r) 1 -"]. 



(1) 



This form is commonly used to describe different re- 
laxation phenomena, including magnetic, optical and me- 
chanical ones, in different complex random systems with 
a distribution of relaxation times i 73 ' 74 However, no gen- 
erally accepted microscopic explanation of this behav- 
ior exists so far. Percolation model puts the restric- 
tions on the possible values of stretching exponents, 
1/3 < 1 — n < IjZS The values of 1 — n below the 
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FIG. 8: Time dependence of the thermoremanent magnetiza- 
tion Mtrm for the LSCNO sample with y=0.25. Mtrm is 
normalized to the value measured immediately after the field 
is set to zero, Mtrm(0). The thick solid line is the best fits 
to a stretched exponential dependence, the dotted thin line - 
to Eq. ([2]) and the dashed one - to Eq. © . 



as pointed out in Ref.[9l[, in Eq. (J2J and Eq. (J3j) it de- 
pends on the time unit used (since this is a fitted value of 
Mtmr at t=i). The dotted line presented in Fig|8]is the 
best fit to Eq. ^ , while the dashed one is the best fit to 
Eq. (|3|). The magnetization is normalized by Mtrm{0), 
the measured value of Mtrm immediately after the field 
is set to zero. The measurement was sufficiently long 
(~ 8 hours) to distinguish between the different possible 
functional forms, given by Eqs. (H])-©. The quality of 
the best fits to the functions J2])-© is not satisfactory, 
as it can be easily seen in Fig. [8] 

The calculated values of % 2 and correlation coefficient 
R 2 in the standard analysis unambiguously show that 
the stretched exponential form ([T]) is the best descrip- 
tion of magnetization decay in LSCNO. The best fit to 
Eq. <[T|) gives 1 - n = 0.32 ± 0.01. The parameter n is 
in perfect agreement with theoretical predictions for the 
SG systems ; 92 ' 93 with the simulations on a 3D Ising SG 81 
and with the experimental results on the canonical SG, 
which also give 1 — nwl/3.— 

IV. DISCUSSION 



mean-field value, 1/3, were reported, albeit they arc be- 
lieved to be non-intrinsio 7 - 6 - or related to fragility of some 
glass formers ^2 The theory predicts that presence of the 
long-range forces (in addition to the obvious short-range 
forces) modifies the glassy relaxation dynamics, reduc- 
ing 1 — n (within the limits given above) j 78 ' 79 The rela- 
tionship ([1]) was theoretically predicted also for the SG 
systemsi^Sr— The numerous experimental data on the 
canonical SG, such as, for example, AgrMn^ 3 - confirm 
validity of this formula for the description of the relax- 
ation of magnetization in these systems^ 9 - The best fit to 
Eq. (P) for LSCN025 is visible in Fig.|S]as the solid line. 

We have also tried to fit the data to other formulas 
used to describe decay of magnetization, a power law 
and a logarithmical dependence, 

M TRM (t) =M t-i, (2) 

M TRM (t) = M - Sln{t), (3) 

respectively^ Eq. has been used to describe relax- 
ation of magnetization in the systems with both long- 
range-AF and long-range-ferromagnetic orderi 84 ' 85 Nu- 
meric simulations predict that this equation should be 
also applicable to the SG systems.— Eq. ([3]) has been 
found to be valid for the systems where the energy bar- 
riers over which magnetic relaxation takes place are uni- 
formly distributed from zero to a certain maximal energy 
and the behavior consistent with this equation has been 
observed in several different SG systems*^— It has been 
noted that when the decay parameter 7 is small, the ex- 
perimental data in the limited time interval can be fitted 
both by Eq. ^ and Eq. ^ equally well.^^ In all these 
equations Mq is one of the fitting parameters but, 



A. Static scaling 

The definitive feature that corroborates the presence 
of the SG phase in LSCNO is the scaling behavior. In- 
troducing the SG susceptibility, XSG, proportional to av- 
eraged (thermally and spatially) 95 ' 96 square of the spin 

correlation function, \sg oc / ^[(ffj ■ Sj) T ]av, allows to 

analyze the SG transition within the framework of the 
second-order phase transition theory, with diverging cor- 
relations at T g and onset of the "order" below T g , and 
next to implement the static scaling hypothesi a 97 ' 98 to 
description of the SGi22 

Experimentally, the Xsg is measurable through 
the dimcnsionlcss nonlinear susceptibility, defined as 
Xnl{T,H) = 1 - M(T,H)/xiH 1 100 - 102 The linear sus- 
ceptibility, Xh comes from the measurements at low 
field. This definition of Xni means that the measured 
temperature-dependent susceptibility, x = M/H, can be 
written as X=Xi(l — Xni)- Note that sometimes the term 
"nonlinear susceptibility" refers in the literature to the 
coefficient at the third power of field in the expansion 
of magnetization in odd powers of field, d 3 M/dH 3 \H=o- 
In the following, by Xni we mean the "whole" nonlin- 
ear part of To use Xni in the scaling analysis, one 
needs to stand up to the problem with estimation of the 
critical region where the scaling should hold and to the 
fact that the state below the SG transition temperature 
is not a state of thermodynamic equilibrium^^ Despite 
these difficulties, Xni has been found to be a good tool to 
investigate a possible thermodynamic phase transition in 
the various 2D and 3D SG systems*!^ 

In FigE] we depicted the susceptibility of LSCN025 
measured at various fields from ^qH=10 Gs to 5 T. The 
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FIG. 9: Magnetic susceptibility of LSCN025 as a function 
of temperature at various fields. The large symbols denote 
ZFC (black spheres) and FC (open squares) data measured 
at fxoH=10 Gs. The solid line is the best fit of CW function 
(plus constant) to the FC data above 12 K. The small open 
symbols represent data measured at 0.02, 0.1, 0.2, 0.5, 1, 2, 3 
and 5 T (from top to bottom) in the FC mode. In the inset, 
the open circles denote the nonlinear susceptibility Xni as a 
function of field at various temperatures at 0.5 K intervals, 
from 6 K to 10 K and from 11 K to 14 K (from top to bottom). 
The solid line is the best fit of the dependence ([5} (see text) to 
the data at T 9 = 10.4 K (denoted as solid circles). The dashed 
lines are the guides to the eye. 



temperature-dependent part of the susceptibility at low 
field and at temperatures above «1.2 Tj is well described 
by the CW dependence, X (T) - Xo = C/{T - 9). The 
best fit to the FC x( T ) data in the range 12 K - 100 
K, showed in Fig|9] as the solid line, gives negative ef- 
fective 6=— 5.3 ± 0.1 K, indicating the presence of the 
AF correlations. We will take only the T-dcpcndcnt part 
of the measured x m to account in the following scaling 
analysis. We assume that x(T) measured at 10 Oe is 
a good approximation of Xi an d thus xi = C/{T — 0). 
At low temperatures, the magnetization increases non- 
linear ly with the field and the xC^O deviates from the 
CW function towards the smaller values at larger fields, 
as it is clearly visible in FigFJ] 

The scaling theory predicts for Xni the relationship 

where t is the reduced temperature, t=(T — T g )/T g , (3 
and 7 are the critical exponents, and /_ (/+) is the scal- 
ing function for t < (t > 0), 101 i 106 i 107 (Another ap- 
proach to the static scaling, sometimes used in the lit- 
erature but clearly giving unsatisfactory results in the 
case of LSCNO, is discussed in Appendix [Bj) The scal- 
ing functions /^(x) behave as const-a; 2 /* in the largc-x 



limit i 100 i 106 This means that right at T g we have 

Xni oc (H 2 ) 1/S , (5) 

where <5 is another critical exponent, related to j3 and 7 
by the scaling law 

(5-1)0=7. (6) 

In the inset to Fig|9] we present the Xni vs H curves 
for various temperatures. They are results of the isother- 
mal cross-cuts of the Xni{T) curves calculated from the 
experimental data for various fields. The position of the 
maximum in the ZFC xCO curve at 10 Oe, equal to 
10.4±0.1 K, is taken as T g . The Xni(H) curves on log- 
log scale change their curvature sign around T g . Right 
at T g , the data are described by a linear dependence on 
the log-log scale and the best fit to Eq. ([5]), shown in the 
inset as the solid line, yields <5=5.8±0.1. 

With this value of 6 we have adjusted j3 and have cal- 
culated 7 from Eq. © to obtain the optimum coincidence 
of the data on two universal curves, one for t < and 
second for t > 0, in the Xnt/\tf vs H 2 /\t\^ plot. The 
best qualitative collapsing of the data to these two sep- 
arate curves has been found for (3=0.75, what implies 
7=3.6 (see Fig[15]in Appendix IB"]) . We have estimated 
the uncertainty of the adjusted parameter to be ^0.05, 
i.e. /3=0.75±0.05 and thus 7 =3.6±0.3. 

To estimate and to visualize the critical temperature 
region where the scaling is valid, it is better to use 
the argument of the scaling function that is linear in 
t. Such improved form of the scaling has been pro- 
posed by Geschwind et al. for the equation Xni = 
H*PfV+i)G(x) = H 2 ^^+^G{H 2 /\tf+^)^. Raising 
the argument x to the power of — 1/(13 + 7) makes it lin- 
ear in t. In the same way, Eq. Q may be reformulated 
as 

The scaling plot with the use of this equation for /3=0.75 
and 7=3.6 is presented in Fig. 1101 In the logarithmic 
scale of y-axis all values of Xn/|^l _ ^j varying over two 
decades, are given equal weight and this allows to notice 
and compare any potential deviations from the universal 
curves at different ordinates. The scaling validity region 
at lower fields is noticeably smaller than at larger fields, 
i.e. the scaling does not work for the same argument 
\t\/H 2 /(P +J ) of function / T for which the scaling at larger 
fields is still valid. In terms of temperature critical re- 
gion, this means that at the field 0.2 T scaling is valid in 
the interval ~ 0.8T g — 1-1T S . The deviations of the data 
outside this T-region from the universal curves are larger 
than the measurements uncertainties (which are of the 
symbol size in Fig. [T0|) . The scaling region expands with 
increasing field and the quality of scaling at 5 T is still 
excellent even at temperatures as far from T g as 0.2T g 
(i.e. the lower limit of our measurement window) and 
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FIG. 10: (Color online) Scaling plot for LSCN025 according 
to Eq. ([7]) with ,3=0.75 and 7=3.6 for the nonlinear suscep- 
tibility xm at various magnetic fields from 0.1 T to 5 T at 
temperatures below (t<0) and above (t>0) the SG-transition 
temperature T 3 = 10.2 K. The error bars are of the symbol size. 



2T g . Such a large scaling region has been found experi- 
mentally to be typical for the canonical SG systems?^ 
in agreement with the numerical simulations for the 3D 
Ising SG, by which the critical region has been estimated 
to extend up to \t\ w 0.554i 

Results of the identical scaling procedure for LSCN017 
and LSCNO50 are similar to those for LSCN025. In par- 
ticular, the scaling according to Eq. (Q| and its improved 
form, Eq. (O, gives clearly the better results than the 
scaling according to Eq. (|B2|) . This remains true even 
for LSCN017 where the strongest AF correlations among 
the investigated samples are expected and thus one might 
expect that the scaling approach described in Appendix 
|B] would work better than that described by Eq. ([7]). The 
best collapsing of the data at different magnetic fields 
onto two branches of the universal curve - below and 
above T g - have been found for LSCN017 for ,3=0.55, 
i.e. the value closer to /3=0.5, predicted by the numeri- 
cal simulations for 3D Ising SG, than to ,9=1, predicted 
by the mean-field theory for isotropic 3D Hciscnbcrg SG. 
The quality of scaling is excellent, as it can be seen in the 
bottom panel of Fig.[TTJ The scaling for LSCNO50 is best 
for /?=0.75, i.e. for the same value as for LSCN025. A 
comparable size of the T-region where the scaling is valid 
starts for LSCNO50 at larger fields than for LSCN017 
and LSCN025, but the quality of the scaling is still very 
good (see the upper panel in Fig. ITT]) . 

The parameter /?, found here as a result of simple ad- 
justing procedure, is a critical exponent for the SG order 
parameter qEA, originally introduced by Edwards and 
Anderson in the model based on the classical (mean-field) 
calculations The corresponding quantum-mechanical 
calculations have been carried out by Fischer^ Since 
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FIG. 11: (Color online) Scaling plot for the nonlinear sus- 
ceptibility xm of LSCN017 (LSCNO50) according to Eq. © 
with ,3=0.55 (/3=0.75) and 7 =3.2 (7=3.8). 



no lattice effects have been taken into account, the ob- 
tained results describe the amorphous SG, where the 
CW law with 6>=0 is found above T SG .^ The possibil- 
ity that the average exchange interaction J is not zero, 
and thus causes that the ferromagnetic order competes 
with the SG phase, has been included into the model 
by Sherrington, Southern and Kirkpatrick (SSK), and 
the formula for extracting qea from the measured x h as 
been given i 96 ! 110 ' 111 The calculations based on the local- 
mean-field approximation suggest that q remains un- 
changed when J changes its signp^i^ The one-component 
SSK model has been modified and extended subsequently 
to describe a two-component Ising-likc magnetic sys- 
tem (with the separate order parameter for each compo- 
nent), where the re-entrant transition to the SG phase, 
both from the ferromagnetic and from the AF phase, is 
predicted ^i^rJJ^ None of these models provides a realistic 
description of LSCNO. Due to lack of analytical expres- 
sion for qEA, the adjusting procedure seems to be the 
best approach to find (3. 
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FIG. 12: The real part of Xac as a function of temperature 
in the vicinity of the transition to the SG phase at various 
frequencies (from top to bottom: 1.00, 2.68, 7.19, 31.3, 138.7 
and 997.3 Hz) for LSCN025. The bottom inset shows the 
peak positions for this sample (open circles) together with 
those for LSCN017 (solid squares) and LSCNO50 (solid dia- 
monds). The thick solid lines are the best fits to Eq. [8] The 
details are described in the text and the obtained parameters 
are given m Table |U In the upper inset, the slowing-down 
formula (Eq. [8} versus Vogel-Fulcher law (Eq. [9j is tested for 
description of experimental data for LSCN025. The thick 
(thin) solid line in left (right) panel is the best fit to Eq. [8] 
(Eq. [9]), correspondingly. Note the logarithmic scale of ab- 
scissa in the left panel and the linear one in the right panel. 



out of equilibrium. According to the dynamic scaling 
hypothesis, the characteristic time r in the vicinity of 
the transition changes with the correlation length £ as 

Fitting to the power law given by Eq.© requires ad- 
justing three parameters: To, T g and the product zv. 
Since T g is the infinitely slow cooling value of Tf (i.e. 
can be regarded as limf^oTf), we carried out a more 
restrictive fit by assigning T g the value of temperature at 
which the ZFC dc-x(T) curve has its maximum. Addi- 
tionally, this allows us to overcome the difficulties caused 
by the small number of the experimental points and the 
relatively large errors compared with the small change of 
Tf within 3 decades of the frequency variation. Having 
T g fixed, it is possible to obtain directly To and zv from 
the linear fits of log r vs log (Tf - T g )/T g . The best fits, 
shown as the thick solid lines in the inset to Fig[l2] (left 
panel) , yield the parameters given in Table HJ 



TABLE I: The dynamic magnetic properties of LSCNO. 



y =0.17 y=0.25 



y=0.50 



T 3 (K) 

TO (s) 
ZV 

E a /k B (K) 
To (K) 
r* (s) 



6.8±0.1 

^q-9.5±0.4 

-6.8±0.4 

-7.6 
6.6±0.1 

^q-6.4±0.2 



10 



10.4±0.1 

-10.5±1.0 



-8.8±1.1 

~ 27 
9.6±0.2 



10" 



S.3±0.4 



13.1±0.1 

^q-9.7±0.5 

-8.7±0.5 

~ 34 
12.2±0.1 

-^q-7.8±0.3 



B. Dynamical scaling 

A more detailed insight into the dynamics of the SG 
state in LSCNO can be obtained from the analysis of 
the Xac measurements. In Figll2j we depicted the tem- 
perature dependence of the real component of Xac for 
LSCN025, measured with 1 Oe ac field amplitude at 
various frequencies. The data have been collected af- 
ter cooling the sample in zero field. The peak position 
Tf moves toward higher temperature with increasing fre- 
quency. This frequency dependence can be described by 
the standard critical slowing down formula, given by the 
dynamic scaling theory,^i& 



To 



T n 



(8) 



In this equation, the characteristic time r describes the 
dynamical fluctuation time scale and corresponds to the 
observation time, t a b s = 1/uj = 1/2tt/, at the tem- 
perature of maximum in x'(T)] to is the shortest time 
available to the system, i.e. the microscopic flipping 
time of the fluctuating entities; Tf is the frequency- 
dependent freezing temperature and v is the critical ex- 
ponent of the spin (or spin-cluster) correlation length £ 
(£ oc [Tf/(T f - T g )] u ). Below T f , the longest relaxation 
time of the system exceeds t b s and thus the system is 



The typical values of tq for the canonical SG, i.e. 
10~ 12 — 1CP 14 s, are of the order of the spin-flip time of 
atomic magnetic moments (^10 -13 s)&Lnrmi As it can 
be scon m Table [fl the LSCN017 exhibits the slowest 
dynamics among the investigated samples. Its charac- 
teristic relaxation time is of the order of ^10~ 10 s and 
is evidently larger than ^10 -13 s expected for the sin- 
gle atomic spins. This strongly suggests the existence of 
spin clusters. Even the shortest To for LSCNO system, 
found for y=0.25, is of the order of ~ 1CP 11 s and thus 
does not exclude the existence of spin clusters, albeit the 
number of spins in the fluctuating entities is expected to 
be smaller, 120 - 121 

The dynamic magnetic properties of a glassy system 
may be tested in the frame of Vogel-Fulcher (VF) lawji2S 



T = T*esxp[E a /k B (Tf-T Q )]. 



(9) 



Taking the LSCN025 data as an example, we show in the 
right panel of the upper inset to Fig JT2] that the change of 
the relaxation time t in LSCNO in the frequency range, 
which has been experimentally accessible to us, can be 
described by this formula equally well as by the power 
law given by Eq.©. 

The problem of discrimination between these two laws, 
the power law given by Eq.© and the VF law given 
by Eq.©, was a subject of debate two decades ago and 
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Eq.© has been found to describe the experimental data 
for the SG systems better than the VF law provided that 
the range of r is sufficiently large.*^ 3 - However, the dif- 
ference is clearly visible only when the variation o f r ap - 
proaches 11 orders of magnitude (see Fig. 2 in Rcf. Il23|). 
This probably explains why the phenomenological VF 
law is still used in literature to describe the frequency 
dependence of Tf in the SG systems. The high-T expan- 
sion of the VF law is identical with that of the power 
law up to terms of order (T /T) 3 . Closer to T (and T g ), 
the VF law can be adjusted to match a power law over a 
large frequency range through the rclatior^ 3 - 



In 



40k B T f 25 
E a zv 



(10) 



With T/-11.6 K found for LSCN025 in the frequency 
range of the experiment, the formula (|10p gives E a /kB ~ 
27 K for this system. Taking this value of E a as granted, 
we have fitted Eq. © to the data. The best fit, marked 
in the right panel of the upper inset to Fig. [12] as the 
solid line, yields r* = l(r 8 - 3±0 - 4 s and T = 9.6 ± 0.2 
K. The analogical analysis for LCN017 and LSCN025 
leads to the parameters given in Table HI The values of 
To are smaller than T g from the power law (Eq. [8]) only 
by a few percent, in accordance with the general trend 
found in the metallic SG systems*^ 

The phenomenological parameter T has been intro- 
duced into the standard Arrhcnius law for glasses to over- 
come the difficulties with keeping the magnitude of r* at 
a physically-meaning level. However, tq from Eq. ([8]) arc 
believed to give more reliable insight into the SG dynam- 
ics than t* does . 120 ' 123 

Despite this, the To can be interpreted as a measure of 
the coupling between the interacting entities.*^ In the 
frame of this picture, To ^ E a /ks indicates a weak 
coupling and T 3> E a /ks a strong one. As it can be 
seen in Table [fl T for LSCNO varies from -0.4 E a /k B 
to —0.9 E a /kB- Thus T is in the intermediate regime 
and confirms the conclusion drawn in Sec. IIII Cl from the 
simple parameter 8 about the presence of some interac- 
tions between the magnetic entities in the system, albeit 
does not settle whether they are single spins or spin clus- 
ters. The energy scale fcsXb is almost equal to E a for 
LSCN017. This betokens the largest coupling between 
the magnetic entities among the examined Ni concentra- 
tions. This is consistent with the observation that \9\ 
in the CW law, describing x(T) in the paramagnetic re- 
gion, achieves its maximal value in the LSCNO system 
for ;y=0.15-0.17 (the inset to Fig. [4]) and suggests that the 
local AF order, partially restored by trapping the mobile 
holes, is strongest at this Ni concentration. 



C. Role of Ni in the Cu-O network 

As it is widely recognized, a 2D dynamic AF order per- 
sists even in the overdoped LSCOji 2 ^ so the dynamical 
AF fluctuations are present in LSC015. They may be 



regarded as a reminiscence of a 3D static AF order in 
the parent compound La2Cu04, destroyed by the frus- 
trating impact of the mobile-hole spins introduced by Sr 
doping*^ Near the famous x=l/8 anomaly, the neutron 
scattering experiments reveal a gap in the low-energy 
spin excitation spectrum, and presence of the magnetic 
incommensurate peaks that has been interpreted as an 
evidence for formation of the static stripes of spins and 
holes in the Cu-0 planes^ 3 -^ In LSC015, this static 
stripe order disappears but the spin gap of 4 meV is 
still visible in the neutron measurements. Introduction 
of even a small number of Ni atoms (y<0.03) into the 
system reduces the spin-excitation energy scaled 

Our study clearly shows that the low-temperature 
phase of Ni-doped LSC015 exhibits all thcrmodynam- 
ical characteristics of the SG phase, including the static 
scaling. The dynamical scaling parameters suggest the 
presence of the ordered magnetic clusters. The XAFS 
measurements reveal that the Ni ions added to the Cu-0 
planes act as the hole absorbers^ 2 - Such Ni 2+ ion with the 
trapped hole has been suggested to form the Zhang-Rice 
doublet^ with the effective spin S'=l/2 and not to dis- 
turb the AF correlations between the S=l/2 spins of the 
Cu 2+ ionsi^ 2 - Our observation of the SG behavior suggests 
that the compensation of the £=1/2 Cu 2+ spin vacancy 
in the Cu-0 network by the effective spin of the complex 
formed by the Ni 2+ ion and the bounded ligand hole is 
not perfect. The observed localized magnetic moment for 
?/<0.09, attributed to the Ni site, is tiny (corresponding 
to 5=0.11) but finite. The T g extrapolated into the SC 
region of the phase diagram takes finite values (see the 
inset to Fig. [3]), what implies that the frustrating effects 
appear in the system with the first added Ni ions. Thus, 
on the contrary to the conclusions of Ref. [42j , our study 
shows that Ni exhibits a magnetic character and affects 
the magnetic correlations between Cu spins in the Cu-0 
planes starting from its smallest concentrations. 

The picture of a glassy system depends on the time 
scale of the used measuring probe. No static component 
in the magnetic signals has been revealed by neutron ex- 
periments on LSCNO up to y=0.029. On the other hand, 
formation of the short-range magnetic order in the vicin- 
ity of the Ni ions is strongly suggested even in the lightly 
doped LSC015 by the zero-field /iSR experiments reveal- 
ing slowing down of the Cu spins fluctuations^^ For 
larger Ni content, the procession of muon spins, indicat- 
ing the long-range magnetic order, has been clearly ob- 
served for y=0.10 at 0.3 K and the static magnetic order 
has been suggested to form even for y>0.07 at temper- 
atures below 2 K, outside our measurement window 
It should be remembered that any relaxation slower than 
10~ 5 -10~ 6 s is seen as a "static" component by /iSRi^^ 
Moreover, the damping rate in the //SR time-spectra 
term representing the muon-spin procession has been 
found to increase with increasing Ni content y for y>0.03 
in the Ni-doped LSCO with a;=0.13 (LSC013). Based on 
this, the less coherent magnetic order, seen by fiSK ex- 
periments in LSC013 with y=0.10, has been suggested 
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to be a precursory state toward the SG stateJ^ Our di- 
rect observation of the SG behavior in Ni-doped LSC015 
is in agreement with this conjecture and our results are 
consistent with the fj.SK measurements. 

As presented in Sec. IIII Bl the susceptibility of LSCNO 
with large y exhibits the CW dependence with the nega- 
tive 9. It is striking that the \9\ reaches its maximal value 
at around j/=0.15-0.17, i.e. where the Ni concentration 
is equal to the hole concentration in the system (see the 
inset to Fig. 2]). This can be qualitatively understood as 
a result of trapping mobile holes by Ni. 

In the framework of this model, trapping the mobile 
holes in LSC015 by Ni ions restores locally the AF order 
in the Cu-0 planesi 42 i 129 The more Ni ions are present in 
the system, the more mobile holes become localized and 
the restored AF order is stronger. This is reflected as the 
increase of absolute value of the effective 9 with y. When 
all mobile holes became localized, i.e. y becomes equal 
to the hole concentration, there is no way to increase the 
AF correlations between the Cu spins and this explains 
why \9\ does not continue to increase with increasing y 
above j/=0.15-0.17. 

As mentioned in Sec. IIII Bl the increase of the calcu- 
lated /i e // per Ni ion with increasing y for y>0.07 means 
that some Ni ions are not in the Ni 2+ L hole-bounded 
state at these concentrations. Thus, a small number of 
the Ni 2+ ions in 5=1 state is probably present in the 
system even for ?/<0.15. However, their influence is not 
able to overwhelm the effect of restoring locally the AF 
correlations between the Cu spins. When the number of 
Ni ions is sufficiently large to localize all mobile holes, 
adding the further Ni atoms to the system means intro- 
ducing the subsequent 5=1 magnetic moments into the 
existing Cu-0 network with the locally AF ordered re- 
gions. This destroys the restored local AF correlations 
between the Cu spins what is observed as decrease of \9\ 
with increasing y above 0.15-0.17. 

The SG transition temperature, T g , in LSCNO in- 
creases linearly with increasing y up to y=0.30, as shown 
in the inset to Fig. |4] In the cluster SG phase of LSCO, 
the monotonic decrease of T g with increasing hole con- 
centration x for 0.03<x<0.05 has been tried to explain 
by the finite-size-scaling hypothesis.— £L According to this 
proposition, T g is expected to be proportional to the size 
of the locally ordered regions L d (d=2), determined by 
the concentration x of the mobile holes destroying the 
AF order, T g oc L 2 oc 1/x, in rough agreement with the 
experimental results 

Applying this model directly to LSCNO would mean 
assuming that the average size of the magnetically or- 
dered region increases with increasing y even up to the 
limit when 1/3 of the Cu ions is replaced by the Ni ions. 
However, decrease of \9\ with increasing y above y~0.17, 
presumably reflecting the destroying impact of the 5=1 
Ni ions on the locally restored AF order, suggests that 
this is not true. Thus, the finite-size effects are proba- 
bly not the only factors determining the value of T g in 
LSCNO up to 7/«0.30. 



At large y, magnetism of LSCNO is presumably domi- 
nated by the Ni ions with the spin 5=1. As it can be seen 
in the inset to Fig. 21 the T g continue to increase with in- 
creasing y above y«0.30, although with the smaller rate 
than for j/<0.30, and exhibits a quasi-linear dependence 
on y up to the largest examined Ni concentration. The 
situation bears some analogy to the standard metallic 
magnetic alloys, where the linear dependence T g on the 
impurity concentration yi m is found at large yi m (in ad- 
dition to the similar T g (xyi m dependence in the regime 
of the interacting single spins at low y im <0.5%, followed 
by the T g cxy^ relationship at larger y)S^ This takes 
place in the interval from about 10 at.% to the magnetic 
percolation limit (where a smooth transition to the inho- 
mogencous long-range order appears). The magnetism 
of the glassy system with yi m within this concentration 
interval is dominated by the large ordered clusters em- 
bedded in the nonmagnetic SG matrix^ In LSCNO, the 
5=1 Ni ions are embedded in the very complex magnetic 
background but probability of forming clusters by these 
ions and the average size of such clusters increases with 
increasing y. Thus, the magnetic behavior of LSCNO 
might be dominated by these clusters at large y. The 
observed increase of the characteristic time for the fluc- 
tuating entities in the system, tq , with increasing y from 
0.25 to 0.50 (see Table [XJ) is consistent with such suppo- 
sition. 

Alternatively, one may interpret the observed SG be- 
havior in LSCNO in the framework of the stripe pinning 
picture. Following this concept, it was suggested that 
the dynamical stripe correlations of spins and holes are 
localized in the vicinity of Ni ions in LSC015 leading to 
the formation of the static stripe order i 127 i 128 

However, the results of the remnant magnetization 
Mtrm measurements, presented in detail in Sec. IIII El 
give us some grounds for speculations contrary to the 
above interpretation. The Mtrm decay in LSCNO 
is described by a stretched exponential function with 
I — ??. « 1/3 exponent, in perfect agreement with the 
theoretical and experimental results for the canonical SG 
systems4ij2Sr^i In "pure" LSNO (without Cu), where the 
presence of stripes has been unambiguously confirmed 
experimentally, the time dependence of the " isothermal" 
remanent magnetization (i.e. obtained after ZFC, apply- 
ing field at 2 K and next switching the field off again), is 
inconsistent with a stretched exponential function.— On 
this ground, it was concluded that the mechanism by 
which the LSNO come back to equilibrium is not the 
same as in the canonical SG, at least in the measured 
LSNO with x=l/3 Sr content where the stripes are com- 
mensurate with the square lattice of the Ni-0 planes, i.e. 
where the charge and magnetic orders have the identical 
periods X 

At the opposite side of the LSCO-LSNO phase dia- 
gram, in the SG phase of "pure" (without Ni) LSCO 
with a;=0.04, the same time dependence of Mtrm as 
that found in LSCNO, and with the same value of 
1 — 77 , w 1/3, has been reported^ However, it has been 
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underscored that the duration of the experiment had not 
been sufficient to exclude the different types of Mtrm if) 
dependence.— In LSCNO, the sufficiently long time of the 
measurement and the standard analysis allow us to ex- 
clude those other possible forms of the Mtrm vs t de- 
pendence (see Sec. MI E]) . 

To recapitulate, the above results can suggest that 
the origin of the SG phase in LSCNO is not related to 
the (potential) presence of the stripes and disorder in 
their array because the mechanism of relaxation seems to 
be identical with the one observed in the canonical SG 
and different from that observed in the reference stripe- 
ordered compound LSNO. 



D. Critical exponents 

In general, the critical exponents are believed to be 
universal and thus capable of describing the behavior of 
the system near phase transition even when the detailed 
microscopic picture is not known. As regards exponent 
zv from the dynamical scaling [Eq.[5], the numerical sim- 
ulations for the 3D Ising SG model give zv = 7.9 (while 
for a conventional phase transition zv = 2 is expected) 41 
Experimentally, the values of zv between 5 and 11 have 
been observed for different SG system s 104 ' 123 and thus 
zv ~ 7 - 9 found for LSCNO (see Table IJ) are inside 
the SG realm and close to the value predicted for the 
Ising-like systems. However, it should be remember that 
it is not possible to distinguish significantly different zv 
values on systems of different spin dimensionality J41 For 
instance, zv ~ 10 — 11 has been reported both for a 
short-range Ising SG system, Feo.5Mn .5Ti0 3 ) 119 i 132 as 
well as for more isotropic vector SG: namely, for a 3D 
XY system, Euo sSri 5 Mn04 f 121 and for 3D Heisenbcrg- 
like AuFe 8% ^3^3. 

With fitted S and adjusted (3, together with the dynam- 
ical critical exponent zv, the other critical exponents can 
be calculated using the scaling equations, in the similar 
manner to 7 from Eq. ([5]). Namely, the heat-capacity 
exponent a can be calculated now from the equation 



a + 2(3 + 7 



(11) 



and the spin-correlation length exponent v can be ob- 
tained from the relationship 

dv = 2-a, (12) 

where d is the dimensionality of the magnetic interac- 
tions. Since the measurements of the canonical SG/metal 
multilayers have showed that even a very weak magnetic 
coupling between layers of 2D SG causes a 3D character 
of the whole system;! 3 ^ we have taken d=3 in this equa- 
tion for LSCNO. The known value of v allows to extract 
the correlation-time exponent, z, from the results of the 
dynamical slowing-down formula fitting. The exponent 
rj that governs the spatial correlation function at T g can 
be deduced from the equation 



The whole set of the obtained critical exponents, in- 
cluding these from the dynamical scaling, is given in Ta- 
ble [IIJ The corresponding values for the 3D Heisenberg 



TABLE II: The critical exponents for LSCNO together with 
the values reported for 3D Heisenberg SG (Ref. [HU) and 
these predicted by numerical simulations for 3D Ising sys- 
tem (first row for a give n exponent - Ref. [8l| , second row - 
Refs.[T3i],[l33] and pi|L 



Heisenberg 



=0.17 



y=0.25 y=0.50 



Ising 



5 


3.3 


6.8±0.1 


5.8±0.1 


6.0±0.1 




0.9 


0.55±0.05 


0.75±0.05 


0.75±0.05 


7 


2.3 


3.2±0.4 


3.6±0.3 


3.8±0.4 


Q 


-2.1 


-2.3±0.5 


-3.1±0.4 


-3.3±0.5 


V 


1.3 


1.4±0.2 


1.7±0.1 


1.8±0.2 


V 


0.4 


-0.23±0.01 


-0.12±0.02 


-0.14±0.02 


zv 

z 


7 

5.4 


6.8±0.4 
4.8±1.1 


8.8±1.1 
5.2±1.0 


8.7±0.5 
5.0±0.7 



6.8 
8-9 
0.5 
0.7-0.8 

2.9 
5.0-6.5 
-1.9 
-4.5- -6.1 

1.3 
2.2-2.7" 
-0.22 
-0.34- -0.40" 
7.9 
6.1 



"Rcfs. 136], 137] and 138]. The rest of static exponents corre- 
sponding to these values of v and rj are obtained using scaling and 
hyperscaling relations. 

and 3D Ising systems are listed for comp arison. The criti- 
cal exponents for AgMn from Ref. |l35 1 , given in the first 
column of Table [TT|. are very typical for the canonical SG, 
being weakly anisotropic Heisenberg-like SG, and thus 
can serve as a reference point for the Heisenbcrg-class 
systems . 121 i 139 i 140 . Since a good experimental realization 
of the Ising SG seems to be difficult,^ - we have used the 
results of the numerical simulations made by Ogiclski,— 



136-138 



as 



(2 - rj)v = 7. 



(13) 



revised by some later large-scale computations 
a reference point for this universality class. The mea- 
sured values for FeMnTiOs, regarded as the best labo- 
ratory realization of the Ising system so far, are roughly 
consistent with these numerical predictions »1^1 

A large variation of the critical-exponent values mea- 
sured for the same SG exists in the literature, par- 
tially because of the different field and temperature 
ranges used^ 9 - Thus any attempt to classify a given 
system to one of the universality classes, based solely 
on these experimental values, should be taken always 
with some caution. In the single crystals of "pure" 
cupratc Lai.9 6 Sr .04CuO4, £=5.9±0.6 and 7=4.3±1.4, 
similar to the values for LSCNO, have been reported.— 
They are significantly larger than those observed in 
the Heisenberg-like SG (see Table [TTJ) and close to 
the exponents measured in the Ising-like systems: 2D 
Rb 2 Cui_ a: Co a: F4 (7=4.5±0.2) and 3D Feo.sMno.sTiOg 
( 7 =4.0±0.3, <5=8.4±0.15). However, the x(T) curve in 
Lai.96Sro.o4Cu04 bifurcates below T g both for H||o6- 
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plane and for H \\c and thus the system has been classified 
ms a 3D Heisenberg-like one4 

Essentially, all exponents lie between those for Ising- 
like and Heisenberg-like SG. The values for y=0.25 are 
practically the same as for y=0.50 and no obvious trend is 
visible with changing y. In the single crystals of "pure" 
nickelatc I^-zSr^NiO^+s (0.02 < x < 0.5) the differ- 
ence between the FC and ZFC magnetization curves has 
been found only for B||a6-plane^ Taking into account 
the aforementioned different result for the " pure" cuprate 
Lai.geSro.o-iCuO^, one might naively expect an evolu- 
tion of the LSCNO properties from the isotropic to more 
anisotropic ones with increasing y. No such clear evolu- 
tion is reflected in the values of exponents for LSCNO. 

Let us note that the critical exponents for LSCNO are 
very similar to those reported for Euo.sSri.sMnO^ clas- 
sified as a 3D XY SG because the bifurcation of x(T) 
curve and frequency dependence of the peak position in 
ac susceptibility has been found only for the field in the ab 
planei^ This suggests that LSCNO may belong to the 
XY SG universality class. However, without measure- 
ments on single crystals or at least on the magnetically 
ordered powders no definite conclusion can be made. 



V. CONCLUSIONS 

The magnetic nature of a Ni dopant in LSCO cannot be 
neglected, even when its concentration is small. The low- 
temperature phase of LSCNO exhibits all the features 
that characterize the spin-glass behavior. The x(^) curve 
displays a ZFC-FC bifurcation and a distinct peak is seen 
in the ZFC branch. A stretched-exponcntial decay of 
the thcrmorcmnant magnetization takes place below the 
transition temperature. The position of a peak in the real 
part of the ac-susceptibility, which is accompanied by a 
step in the imaginary component, is frequency-dependent 
and this dependence is described by the standard critical 
slowing-down formula. The characteristic time that gov- 
erns the internal dynamics of the system indicates pos- 
sible existence of clusters. The nonlinear part of the dc 
susceptibility exhibits scaling behavior characteristic for 
the SG. All these features - when taken together - confirm 
presence of the SG state in LSCNO at low temperatures. 
The critical exponents lie between those characteristic for 
Ising-like and Heisenberg-like systems. 

The transition temperature decreases linearly with de- 
creasing Ni content y and extrapolates to K at y=0, 
suggesting that Ni exhibits a magnetic character and af- 
fects the AF correlations between Cu spins in the Cu-0 
layers starting from the smallest concentration y, in the 
superconducting region of LSCNO phase diagram. It is 
possible that the (cluster) SG phase coexists with the lo- 
cal AF order partially restored by trapping holes at Ni 
sites. 
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FIG. 13: The average distance between Cu/Ni and apical 
oxygen atom in LSCNO with various Ni content, y. The 
solid circles mark the direct results of Rietveld refinement 
procedure while the open circles denote the same results with 
the assumption that zo(2) = 0.181 A= const. The solid line is 
the linear interpolation of the results marked as open circles. 
The small open triangles mark the data from Ref. 14511. th e 
large solid diamond - the value of -Rcu-0(2) from Ref. |142( |. 
and the open star - i?Ni-0(2)(y= 0.06) from Ref. [ij. The 
arrow indicates -Rcu/Ni-0(2) (y= 0.06), marked as the solid star 
and calculated as described in the text. Inset: the refined z 
coordinate of 0(2) apical atom. The solid line marks the value 
0.181. 
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Appendix A: Apical oxygen atoms 

The Rietveld refinement procedure allows to determine 
the values of free atomic coordinates. The position of 
the apical oxygen atom 0(2) in the tetragonal unit cell 
of LSCNO is (0,0,z). The zo(2) values as a function of 
Ni content are shown in the inset to Fig. [T3] The X- 
ray diffraction measurements arc not light-atoms sensi- 
tive; therefore the uncertainties of zo(2) are relative large. 
Within the error margins, there are no obvious trend in 
the zo(2) vs V dependence and zq(2) can be accepted as 
being constant and equal to 0.181(3) in the whole in- 
vestigated doping range. The uncertainties in the deter- 
mined zq(2) are transferred onto the calculated average 
interatomic Cu/Ni-0(2) distances, i?cu/Ni-0(2) = Ro(2)- 
Thus, the found i?o(2) vs V behavior mimics zo(2) vs V 
dependence, as it can be clearly seen in the main panel 
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of Fig. [13] However, taking zo(2)=const as granted re- 
sults in smooth linear dependence of i?o(2) vs V (see 
open circles in Fig. [T3|) that mimics linear decrease of 
c with increasing y. This means that within the error 
margins our measurements do not reveal a nonmono- 
tonic i?o(2) vs y behavior observed in Ref. (4^|. As 
mentioned in Sec. IIII A[ XPD is not atom-specific in the 
sense that only averaged Cu/Ni-0(2) distance is mea- 
sured. For a given y, this measured i?o(2) distance can 
be expressed as y ■ Rm~o(2) + 0- - y) • #Cu-0(2)> where 
Rm-0(2) (i*Cu-0(2)) is Ni-0(2) [Cu-0(2)] distance, re- 
spectively. Ro(2) for y=0.06 can be calculated in this 
way with the use of J?Ni-0(2) =2.250(12) A obtained for 
this Ni content by Haskel et al. from EXAFS spec- 
tra (Ref. H3) and i? Cu _ 0(2) =2.397(5) A from our XPD 
measurements for undopcd (y=0) structure. The calcu- 
lated i?o(2)(y=0-06)=2.392 A is in perfect agreement 
with the value 2.391 A obtained directly from Rietveld 
analysis with the assumption that zo(2) = 0.181 A (sec 
Fig. [T3")) . No significant differences were found in the 
^Ni-0(2) values for y < 0.07. 42 ' 47 The linear dependence 
of i?cu/Ni-o(2) vs y, when extrapolated to j/=l, yields the 
value that is 0.04 A larger than i?Ni-O(2)(z/=0-06) (see 
the solid line in Fig. [T3|) . This fact gives some grounds 
for speculation that i?Ni-o(2) mav change for y > 0.07. 



Appendix B: Scaling analysis 

Since the different approaches to the static scaling in 
SG exist in literature^ ' 103 ' 108 it is worth to explore them 
one after the other to check, which of them describes 
the experimental data in the best way. In the follow- 
ing, we will present the details of scaling analysis for 
LSCN025, but we have also carried out the similar anal- 
ysis for LSCN017 and LSCNO50. 

In the derivation of the Eq. 2] the interactions in the 
system are assumed to average perfectly to zero. 95 ' 101 ' 106 
Since evidently this is not our case and 9 in the CW 
law describing xi °f LSCNO has a substantial value, 
we first tried to follow the approach based on expand- 
ing magnetization M in odd powers of xiH (instead of 

/j ^102.103.143.144 

M(T, H) = xiH - a 3 ( X iH) 3 + a 5 ( X iH) 5 - a 7 ( X iH) 7 + 

(Bl) 

The coefficients 03, as,... here are T-dependent. Since 
now xi has not a simple Curie form, xi ex ^-/T, the 
scaling described by Eq. Q is replaced by the following 
relationship 1 ^ 



Xnl 



9t 



XiH 
\t\ 2 



(B2) 



where g- and g+ is another pair of the scaling functions. 

Since both pair of scaling functions, f T (x) from Eq. Q 
and g^(x) from Eq. (|B2[) . behave as const-cc 2 / 5 in the 
large-a; limit ) 100 ' 103 ' 106 the same value <5=5.8±0.1 found 
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FIG. 14: (Color online) Scaling plot for LSCN025 according 
to Eq. (|B2[) with /3=0.75 and 7=3.6 for the nonlinear suscep- 
tibility Xnl at various magnetic fields from 0.1 T to 5 T at 
temperatures below (t<0) and above (t>0) T 3 =10.2 K. 



from the power dependence of xm vs H at T g (Eq. (J5J); 
see also the inset to Fig. ^ is used in both scalings to 
calculate 7 from the adjusted /3 [Eq. |JBJ]. To perform 
scaling according to Eq. (|B2[) , /? is adjusted in such a 
way that coincidence of the data on two curves in Xni/\t\ 
vs xiH /\t\^ +1 ^ 2 plot is as good qualitatively as possi- 
ble. The best results have been obtained for /3=0.75 and 
7=3.6. As it is seen in FigfTJl the quality of the scaling 
is not satisfactory. 
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FIG. 15: (Color online) Scaling plot for the nonlinear suscep- 
tibility Xnl of LSCN025 according to Eq. @ with /3=0.75 and 
7=3.6 and at the same fields and temperatures as in Fig. 1141 
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To improve the quality of scaling we have employed 
Eq. (U]). The best collapse of the data onto the two uni- 
versal curves has been obtained once again for /3=0.75 
with the uncertainty estimated to be ~0.05. This implies 
7=3.6±0.3. The resulting scaling is shown in FigfTS] The 
quality of scaling is much better than this obtained by 
employing Eq. (|B2|) and almost perfect for fields larger 
than 0.2 T. Thus, this form of scaling [i.e. given by 



Eq. g])] is used for LSCNO. 

However, while the log-log plots presented in Figs. [14] 
and [15] are enough to make direct comparison between 
the quality of obtained the best scalings according to 
Eq. (|B2[) or Eq. ((4]), to estimate the critical temperature 
region where the scaling is valid, it is better to use Eq. (0 , 
in which the argument of the scaling function is linear in 
t, as described in Sec. IIV Al 
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